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Abstract — Time-delay vibration Control Technology is 
applied to the precision instrument workbench, whose influence 
is analyzed jointly by the control of fractional derivative and 
delayed feedback. Delay feedback control with fractional 
derivative is considered the workbench modeling. The stability 
switch theory is used to analyze the stability of the delayed 
dynamical systems, and find the system stable region. In the 
frequency domain, the objective optimization function is 
established based on the frequency domain. The delay feedback 
control parameters are optimized by the PSO (particle swarm 
optimization) and we can get the optimized control parameters. 
What’s more, the system vibration response simulation analysis 
is taken under time domain. Simulation results show that the 
fractional derivative and delay vibration control technology 
combined with the active suspension system can achieve good 
damping effect under harmonic excitation and random road 
excitation. Simulation results show that fractional derivative 
and delay vibration control technology combined with the 
workbench is less than the without fractional derivative and 
delay vibration control technology combined value is reduced by 
99%, The value of the velocity and acceleration feedback is 
reduced by 99% and 68% under harmonic excitation. It shows 
that the fractional derivative and delay vibration control 
technology combined with the precision instrument workbench 
can effectively reduce vibration and provide a theoretical basis 
for active control technology. 

Index Terms — Fractional; Workbench System; Time-delay; 
Stability; Particle Swarm Optimization (PSO) 

I. Introduction 

In the aerospace, automotive industry, instrumentation, 
weapons, construction and engineering machinery and other 
fields, precision instrument workbench has a wide range of 
use. As the external environment and equipment work always 
produces vibration, vibration on the performance of the 
instrument causes bad impacts, so the normal use of the 
equipment and vibration-damping technology which is 
closely related to effective ways to reduce vibration is 
particularly important. 

In recent years, the theory and application of time-delay 
control and fractional-order calculus have become a hot topic 
at home and abroad. They have attracted the attention of many 
research institutions and scholars, and have obtained a large 
amount of research results. In 1992, Olgac, McFarland and 
Holm-Hansen [1] first proposed the concept of time delay 
dynamic vibration absorber. Subsequently, Olgac [2] has 
designed a dual-frequency time delay dynamic vibration 
absorber which can be given a completely fixed 
dual-frequency inhibiting the vibration of the main system. 
Jalili and Olgac [3] solved the stability of the time-delay 


vibration system by using the method of graphic setting. 
Olgac , Renzulli [4] and many people realize the time delay 
dynamic vibration absorber can achieve parameter real-time 
control. Xu and Zhao [5-6] studied the time delay nonlinear 
dynamic vibration absorber and the nonlinear dynamic 
vibration absorber. The result shows that the stability control 
range of the nonlinear time delays, dynamic vibration 
absorber is better than that of the linear time delay dynamic 
vibration absorber wider range. Zhang Wenfeng, Hu Haiyan 
[7] and others studied the time delay factors in the vehicle 
suspension 'ceiling' damping control. Compared with the 
integer order, fractional calculus is the main advantage in 
describing the system memory and heredity. Deng and Li [8] 
studied the robust stability of fractional systems with multiple 
delays. Mihailo and Lazarevic et al. [9-10] used 
Bellman-Gronwall theorem to solve the stability problem of 
time-delay fractional-order system in a finite time interval. 
Chen and Moore [11-12] used the Lambert W function to 
analyze the stability of fractional order systems with time 
delay. Based on the principle of the argument and the Hassard 
principle of the time-delay system, the stability criteria of the 
fractional-order dynamic system with time delay are proposed 
by Wang et al. [13] This paper provides a theoretical basis for 
the application of fractional derivative and active damping 
control technology to the damping of the precision instrument 
table. 

In this paper, we study the influence of fractional derivative 
and time delay damping control on the vibration of a precision 
instrument table with active damping control technology. The 
purpose of this paper is to improve the performance of the 
workbench more comprehensively by increasing the number 
of control parameters, adding time-delay state feedback with 
fractional derivative in the workbench system, establishing 
frequency-domain characteristics of the ground excitation 
according to the table acceleration. The objective function of 
the system vibration is obtained. The control parameters are 
selected through optimization analysis. The optimal feedback 
coefficient, the time delay and the fractional order are 
obtained, and the simulation is carried out in the time domain. 

II. Mechanical model 

Considering the structure of the workbench system and the 
form of vibration, neglecting the degrees of freedom in the 
pitch and roll directions of the system, only consider the 
vertical vibration workbench. The system is simplified as a 
two-degree-of-freedom system consisting of a workbench and 
a vibration absorber. 
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Tab. 1 workbench model parameters 


quality 
of the 
platform 

quality 
of the 
vibration 

absorber 

stiffness co 
efficient 
of the 
platform 

stiffness coeffi 
cient of the 
vibration 

absorber 

damping 
coefficient 
of the 
platform 

damping 
coefficient of 
the vibration 

absorber 

m u kg 

m 2/ kg 

k\/N. in' 

ki/N. rri 1 

C\/N. s. in' 

Ci/N. s. in' 

150 

15 

15000 

1500 

30 

3 


According to the purpose of the study, simplify the 
workbench, and establish the mechanical model as shown 
below: 



Fig.l. 1/2 workbench model with time delayed feedback 
control based on fractional derivative 

Let x 1? x 2 be the generalized coordinates, and obtained the 
differential equations of the system. 

m,!, + 2q (i, - q) + 2k } (jq -q)-k 2 (x, - jq) - c 2 D a (x 2 - jq) - gjq (t- r) = 0 

< ( 1 ) 
m 2 x 2 + k 2 (x 2 - jq) + c 2 D a (x 2 - jq) + gjq (t- r) = 0 

ml is the quality of the workbench, and m2 is the quality of 
the vibration absorber. k 1? ci, respectively, is the 
stiffness coefficient and damping coefficient of the 
workbench. k 2 , c 2 , respectively, is the stiffness coefficient and 
damping coefficient of the vibration absorber, g is the gain 
coefficient of displacement feedback, x is delay. D a is the 
fractional differential operator. The workbench specific 
parameter values as showed in figure 1. The design of 
dynamic vibration absorber parameters is based on the basic 
principle of the dynamic vibration absorber. When the 
dynamic vibration absorber resonates with the external 
excitation frequency acting on the main system, the dynamic 
absorber absorbs energy to make the main system have the 
best vibration-damping effect. Assume that the systems are 
motivated by the periodic excitation, ql= sin (14t). x 1? x 2 , 
respectively, is the vertical displacement of the workbench 
and the vibration absorber. 

III. System stability analysis 
A. System stability conditions 

The Laplace transform is applied to the differential system 
equations: 


Equation (3) is transformed into as follows: 

CE(s 1,n ) = P(s 1,n ) + Q(s 1 "' ) e“ rs = 0 (4) 

n is the common denominator of the exponential power of 
complex variables. The function of P, Q as follows: 

P{ s l,n ) = m 1 m 2 x 4 + 2m l c l s 3 + (2 m 2 k l + m 2 k 2 + m { k 2 )s 2 + 2 c x k 2 s 

+2k x k 2 + (c 2 m 2 + m 1 c 2 )v + 2 c x c 2 s + 2 k x c 2 s 

Q(s 1/n ) = -m 2 gs 2 + gk 2 

According to the system stability and sufficient conditions 
when all the characteristic roots s t (i = 1,2,..., n) of the system 
have negative real parts (in the left half plane of the [s] plane). 
That is \arg(Si)\> tt/ 2. The system is stable, as showed in figure 
2 . 



Introduce an auxiliary complex variable 2-S\!n, and the 
characteristic equation with the fractional power in s domain 
is transformed into the characteristic equation of integrative 
power in 2 domain. The characteristic equation is written as 
follows: 

CE{A) = P(A) + Q(A) e“ a " = 0 (5) 

According to the characteristic equation in k complex plane, 
the system can obtain the stable region, \arg{s l )\>n!2. As 
showed in figure 3. 



It can be seen that the stability region of the fractional-order 
system is larger than the stability region of the integer-order 
system. Therefore, the system control parameters have more 
intervals for adjustment. The system has better robustness. 



X 


2 c { sQ +2 k x Q 


X 


0 


( 2 ) 


m/ + 2c x s + 2 k { +k 2 +c 2 s a - ge~ TS -k 2 - c 2 s a 

-k 2 - c/ +ge~ TS m 2 s 2 + k 2 + c 2 s a 

Where, X 1? X 2 , and Q, respectively, is the Laplace 
transforms of x 1? x 2 , and q. 

To get the characteristic equation of the system: 


CE( s) = 


niyS 2 + 2c x s + 2k x + k 2 + c 2 s a - ge 


TS 


a 


-ts 


k 2 c 2 s 
m 2 s 2 + k 2 + c 2 s 


-0 


( 3 ) 


-k 2 - c 2 s +ge~ 

Including, ae (0, 1), because of continuing e‘ TS beyond the 
exponential function, analyzed the stability of the system 
using stability switch method. 


B. System Stability Criteria 

Firstly, the Rekasius substitution is used to transform the 
characteristic equation containing the transcendent term into 
an ordinary algebraic equation, the Rekasius as follows: 


_ l-Ts _ l-TA n 
e ~ 1 + Ts ~ 1 + TA n 
2 , 

r — — (tan ( coT ) + k7i ) k e □ , T e/? 

CO 


( 6 ) 

( 7 ) 


By equation (5) and (6), we can obtain the characteristic 
equation as follows: 

Ce { X ) = P(A)(1 + TA n ) + Q(^L)(l - T A n ) (8) 

Let s=jco, then 2=co lln Q JJT/271 . Substitute 2=co 1/ne j n/2n into the 
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equation ( 8 ), and separated the real part and imaginary part. 
And get the equation. 

Ce(co l/n e j7r/2n ) = Ce m (■ co 1/n e j7r/2n ) + jCe 3 (■ co l/n e j7r/2n ) (9) 

Let the real and imaginary parts equal zero. We can get the 
simultaneous equations of T and co. 


Ce^(co l,n e j7r/2n ) = 0 

,1/n „jn/2n ^ — Q 


( 10 ) 

Ce 3 (a 1,n e’~""' ~ 

If there is not a real solution of equations (10), the delay of 
the motion trajectory of change will not make the system 
characteristic root to the imaginary axis; if there is a real 
solution, existing system characteristic root is located in the 
imaginary axis. 

By equations (7) all the critical time delay system can be 
obtained. 

The traversal direction of the root trajectory of the system 
at the critical characteristic root is judged by the characteristic 


root movement trend RT 


s=joo\-> 


the characteristic root 


movement trend is denoted as follows: 


RT 



f 

f 


\ 



ds 



. - sgn 





s=jeo i ° 


dr 

s=ja>i 



V 

V 



( 1 ) 


Where, 


ds 


dr 


sQ(s)e 


S = JCOi 


P'(s) + (Q'(s)-TQ(s))e 


s=ja>i 
T = T: 


If RT | s=jcOi=- 1, the number of the system is reduced by 
two unstable roots. If RT \ s=]co{=+\, the number of the system 
is increased by two unstable roots. 

When the delay is zero, using the equation larg(si)l> 7 i /2 to 
judge whether the system is stable. At the time delay is not 
equal to zero. System stability is closely related to the solution 
of equations (10). If equation (10) does not have a real 
solution, the system is stable with full time delay. If the system 
( 10 ) has a real root, the system will have a critical 
characteristic Sf=jcOi on the imaginary axis, and there exists a 
critical time lag accordingly z ik = (2/co,) (tan 1 (cofT) +kn), and 
then at the critical delay, according to the trend of the 
characteristic root movement RT\ s = j(0i judge whether the 
system trajectory through the imaginary axis, so as to judge 
the critical time delay interval of the system stability. 

IV. Control Parameter Optimization 

Particle Swarm Optimization (PSO) [14-16] is an adaptive 
stochastic optimization algorithm based on search strategy. 
As the optimization algorithm is simple and effective, the 
convergence speed is quick, and the parameter setting is few. 
There is no centralized control constraint, and it has strong 
robustness of the system. So it is widely used in 
multi-objective function optimization and intelligent system 
control. The particle swarm optimization (PSO) algorithm is 
combined with the time-delay fractional order damping, and 
the amplitude-frequency characteristic function of the 
workbench system is used as the objective function, and 
solves the optimal value. 

According to the needs of the actual algorithm, the 
three-dimensional search space is established to represent the 
gain, time delay and fractional order respectively. The 


particle swarm is generated separately, considering the large 
gap between the parameters to be optimized X= (XI, X2,*** 
Xn), Y=( Yi,Y 2 ,...Y w )> Z=( Zi,Z 2 ,*”Z n ). The particles of the 
search space are denoted as vectors, respectively, X,-= (x ih x i2 ) 
, Y/= (y, 7 , y i2 ) , Z i= (z/ 7 , z i2 ) , and calculate the fitness value 
corresponding to each particle position. The velocity of the 
i -th particle is V, = (V, 7 , \ i2 ) . During the each iteration, the 
particle updates its speed and position through the extremity 
of the individual and the global extreme, updating the 
formula: 


v£ = ®v*-‘ + 


eft ( pbest id -x k id 1 ) + c 2 r 2 ( gbest id - x* 1 ) 


y K — Y 

^ id id 


k -1 


+v 


k -1 


id 


Where, co is inertia weight, cl and c2 are acceleration 
constants. Adjust the maximum learning step, rl and r2 are 
random functions to increase the randomness of a random 
search. In the process of optimization, in order to avoid the 
blindness of searching the optimal particles, according to the 
optimum conditions and actual conditions of the objective 
function, set the reasonable range of gain, delay and fractional 
order. First use maple to calculate IH (co) I, and then bring into 
Matlab using particle swarm optimization algorithm and get 
feedback gain g and delay x and a. 

In order to facilitate the numerical analysis, and introduced 
the dimensionless parameter. 


,A r * = JA«,=^,/?=Ar = k c= ^ g s ai= i ,v i 

rrur M 0 Co pk x m 1 k x \l m 


m x t 


( 12 ) 


The differential equation is written as a Laplace transform. 

X/ + 2cs(X, - Q) + 2(X, - Q) - y(X, - X,) + 0cs° - (X, - X,) - " = 0 

(0 

X,s 2 + Pa,c-(X, - X, )s° - ya, (X, - X,) + g,Xf‘ = 0 

(O' 

Where, .v=jco. 

The optimization objective function is Amplitude - 
Frequency characteristics of workbench for pavement 
excitation. 


min 


H(co) 


vH 


(2 + 2cs) 

V / 


Q(oj) 


Ai4 — 42^1 


( 2 ) 


s.t< 


g>0 

0 < T < 1 
0 < a < 1 


Where 


Ai~ s + ^ jCS + 2 + y — fics a — g k e TS , A n - fie — s a - y 

co 


4i =ya l +g k e~ TS -Pa x c — s 

co 


a 


1 

A 99 = pay -s a - ya x + V 
co 


The result is shown as follows: 

Optimal individual fitness 



Fig. 4 Fitness function iterative optimization 
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By analyzing the impact on the objective function, we can 
get a set of optimization parameters a = 0.4, g = 9986.1 N/m, 
T = O.U. 

Providing that the time delay is zero, the characteristic 
equation of the fractional time-delay system degenerated to 
the characteristic equation of the general fractional-order 
system. To make the system stable, it should be met 

\arg(l i )\>n/2n=n/4 (14) 

Where n=a l = 2.5. 

Substituting a = 0.4, g = 9986.1 N / m, r = 0.1s into the 
characteristic equation (5) and get 

\ 2 =6.34151/5.1125 > \ 4 =1.41821/2.4182 
A 56 =-6.01521/6.1525 > = -2.2319±/1.1255 

It is proved that the characteristic root satisfies (14). The 
system is stable when the time delay is zero. 

When the time delay is 0.1s, the optimized control 
parameter a - 0.4, g - 9986.1 N / m is substituted into the 
system of equations (10). Use the Maple software to calculate 
the solution (T, co ), where the value of co as follows: 

co l2 =-5.3414 ±/2.2518 c% A = -69.8681 ±/15.6601 

co 56 =7.4253 ±i \.2060 co ls =69.2512+ Z5.5239 

f 

We can see the real domain without solution in equations 
(10), the system does not have the characteristic root of 
critical stability. Correspondingly, there is no stable critical 
time delay. This indicates that the change of time delay does 
not cause the trajectory to traverse the imaginary axis of the 
complex plane, and does not cause the system stability to 
change. Since the system is stable at zero time delay, the 
system is stable over the entire time-delay interval. This 
shows that the selected optimal control parameters are 
feasible. 

V. Vibration response analyses 

In order to test the effect of fractional order delay damping 
in the system vibration reduction, and by excluding the 
workbench of a class of linear fractional order vibration 
absorbing system for comparison. Analysis of frequency - 
domain characteristics of time - delayed fractional - order 
absorber system and time domain response of vibration 
damping effect. 

A. Frequency Domain Characteristics 


To optimize parameters of a = 0.4, g = 9986.1 N / m, r = 
0.1s into the amplitude, frequency characteristics of vibration 
response function, and get the figure 5. 



From the figure, in the frequency range of 4 - 15Hz, 
because of the time delay, it has a damping effect. It can 
change the vibration phase, thereby reducing the vibration 


amplitude of the main system. Fractional state feedback can 
not only adjust the damping force, but also adjust the elastic 
restoring force. At the same time, it can adjust the control 
force size and phase, so that the system to achieve the optimal 
damping effect. The general advantages of active vibration 
absorber with fractional time delay dynamic absorber. That 
can achieve real-time tuning. When the passive dynamic 
vibration absorber reduced vibration absorption effect due to 
the system frequency changes, the fractional time delay 
dynamic vibration absorber can be tuned, so that the main 
system to achieve the minimum vibration. From this can be 
drawn that compared with the passive system, the active 
system with fractional order time-delay feedback control can 
reduce the vibration of the workbench. If we reasonably 
control the delay and fractional size, it will match with the 
feedback system. The vibration of the system can be 
significantly reduced. It is shown that the time-delay 
fractional-order dynamic vibration absorber has a great 
advantage over the traditional passive vibration absorber in a 
frequency band. 

B. Time Domain Simulation 

Passive feedback systems without fractional order and time 
delay are respectively compared with the active feedback 
system with fractional order and time delay. Analyze the time 
domain simulation of vibration response under the action of 
harmonic excitation (q=sin (14t)). The results are shown in 
Figure 6. 



a) Workbench displacement response without time delay 

and fractional 



b) Workbench displacement response with time delay and 

fractional 



time (I) 

c) Workbench velocity response without time delay and 

fractional 
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d) Workbench velocity response with time delay 

delay and fractional 



e) Workbench acceleration response without time delay and 

fractional 



f) Workbench acceleration response with time delay and 

fractional 

Fig.6 The simulation comparison chart of the workbench 

stability index 


The simulation shows that in the frequency range of 4 ~ 
15Hz, the time-delay fractional-order dynamic vibration 
absorber has a great advantage over the traditional passive 
vibration absorber. Including delay fractional order feedback 
displacement decrease 99% than without delay and fractional 
order feedback, the vibration of the main system is almost 
completely absorbed, the velocity feedback is reduced by 
99% and the acceleration feedback is reduced by 68%, as 
showed in table 2. This indicates that when the system is 
stimulated by the external incentive, the fractional-order time 
delay dynamic vibration absorber has the ability to make the 
external excitation fast decay performance, and the 
attenuation effect is very obvious. Therefore, the fractional 
derivative and delay damping active control technology is 
applied to the damping of the precision instrument workbench, 
controlling the feedback gain, delay, differential operator 
value, can effectively reduce the vibration. 


Table2. The workbench stability index 


Vibration response 
values 

displacement 

feedback 

(mm) 

velocity 

feedback 

(m/s) 

acceleratio 

n feedback 
(m/s' 2 ) 

Without delay and 
fractional 

feedback 

0.7525 

7.8579 

0.0536 

With delay and 
fractional 

feedback 

0.0023 

0.0022 

0.0174 


VI. Conclusion 

In this paper, the theory of time-delay fractional-order 
dynamic vibration absorber is used to construct the model of 
precision instrument workbench. Fusing the theory of 
fractional - order dynamic vibration absorber and particle 
swarm optimization. The time - delay and fractional - order 
feedback control is optimized and the model is simulated 
numerically. The following conclusions are obtained: 

Firstly, for the precision instrument workbench system, the 
active feedback control is introduced by using the time - delay 
and fractional - order dynamic vibration absorber theory. The 
stability model of the workbench is established and stability 
analysis is done by using the stability switching principle, and 
obtained the system stability interval, and the results show 
that the optimization results satisfy the system stability 
condition. 

Secondly, based on particle swarm optimization (PSO) 
algorithm, the time - delay and fractional - order vibration 
parameters are optimized to optimize the objective function 
with the amplitude - frequency response of the workbench. 
The simulation shows that active systems with 
fractional-order time-delay feedback control can be tuned in 
real time, and adjusts the phase and control the size. Under the 
simple harmonic excitation, in the frequency range of 4 ~ 
15Hz, including delay fractional order feedback displacement 
decrease 99% than without delay and fractional order 
feedback, the vibration of the main system is almost 
completely absorbed, the velocity feedback is reduced by 
99% and the acceleration feedback is reduced by 68%. It is 
shown that the time-delay fractional-order dynamic vibration 
absorber has a great advantage over the traditional passive 
vibration absorber. It can significantly improve the vibration 
absorber damping performance, so that the main system can 
achieve the minimum vibration. 

Third, it is a relatively new research method to apply the 
fractional calculus theory and time delay damping technology 
in the vibration reduction project of the precision instrument 
table, and it is also the main innovation of this paper. The 
results of this paper prove that the method is feasible and 
effective in theory, and has certain engineering value. 
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